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I. INTRODUCTION 



Polarized M0ller scattering has been a well-studied low-energy reaction for close to eight 
decades now but has attracted especially active interest from both experimental and theo- 



retical communities due to the recent rapid progress in measuring spin-dependent observables. 
Since the nineties the interaction has allowed the high-precision determination of the electron- 
beam polarization at SLC Q, SLAC Q (j, JLab |jj and MIT-Bates Q. A M0ller polarimeter 
may also be useful in future experiments planned at the ILC [tJ. In addition, polarized M0ller 
scattering can be an excellent tool for measuring parity- violating (PV) weak interaction asym- 
metries (8j. 

The first observation of Parity Violation in M0ller scattering was made by the E-158 ex- 
periment at SLAC [||, which studied M0ller scattering of 45- to 48-GeV polarized electrons 
on the unpolarized electrons in a hydrogen target. Its result at low Q 2 = 0.026 GeV 2 , 



Ap V = (1.31 ± 0.14 (stat.) ± 0.10 (syst.)) x 10 7 [lOj allowed one of the most important 



parameters in the Standard Model - the sine of the Weinberg angle - to be determined with an 
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accuracy of 0.5% (sin 2 $w = 0.2403 ± 0.0013 in the MS scheme). A very promising experiment 
measuring the e-p scattering asymmetry currently running at Jefferson Lab, Qweak jllj ]. aims to 
determine sin 2 9w with relative precision of 0.3%. The next-generation experiment to study e-e 
scattering - MOLLER, planned at JLab following the 11 GeV upgrade - will offer a new level 
of sensitivity and measure the parity-violating asymmetry in the scattering of longitudinally 
polarized electrons off an unpolarized target to a precision of 0.73 ppb. That would allow a de- 
termination of the weak mixing angle with an uncertainty of ±0.00026 (stat.) ± 0.00013 (syst.) 
fl2| . or about 0.1%, an improvement of a factor of five in fractional precision when compared 
with the E-158 measurement. 

Since M0ller scattering is a very clean process with a well-known initial energy and low back- 
grounds, any inconsistency with the Standard Model will signal new physics. M0ller scattering 
experiments can provide indirect access to physics at multi-TeV scales and play an important 
complementary role to the LHC research program [l3| . 

Obviously, before we can extract reliable information from the experimental data, it is nec- 
essary to take into account higher order effects of electroweak theory, i.e. electroweak radiative 
corrections (EWC). The inclusion of EWC is an indispensable part of any modern experiment, 
but will be of the paramount importance for the ultra-precise measurement of the weak mixing 
angle via 11 GeV M0ller scattering planned at JLab. In general, from the theory point of view, 
the interpretability of e-e scattering is exceptionally good. However, to match the precision of 
MOLLER experiment, theoretical predictions for the PV e-e scattering asymmetry must in- 
clude not only full treatment of one-loop radiative corrections (NLO) but also leading two-loop 
corrections (NNLO). 

A significant theoretical effort has been dedicated to one-loop radiative corrections already. 
A short review of the literature to date on that topic is done in [lij]. In [lH, we also calculated 
a full gauge-invariant set of the one-loop EWC both numerically with no simplifications using 
FeynArts [IBj], FormCalc [lfij], Loop Tools [l6| and Form [l7[ as the base languages and by 
hand in a compact form analytically free from nonphysical parameters. The total correction 
was found to be close to —70%, and we found no significant theoretical uncertainty coming 
from the largest possible source, the hadronic contributions to the vacuum polarization. The 
dependence on other uncertain input parameters, like the mass of the Higgs boson, was below 
0.1%. 

It is possible that a much larger theoretical uncertainty in the prediction for the asymmetry 
may come from two-loop corrections. Paper [l8[ argued that the higher order corrections are 
suppressed by a factor of either about 0.1% or 5% (depending on a type of corrections) relative 
to the one-loop result. However, since the one- loop weak corrections for M0ller scattering are so 
large and since the 11 GeV MOLLER experiment is striving for such unprecedented precision, 
we believe it is now worth looking into evaluating two-loop weak corrections. 

One way to find some indication of the size of higher-order contributions is to compare results 



that are expressed in terms of quantities related to different renormalization schemes. In [19 



we provided a tuned comparison between the result obtained with different renormalization 
conditions, first within one scheme then between two schemes. Our calculations in the on-shell 
and CDR schemes show a difference of about 11%, which is comparable with the difference 



of 10% between MS (20[ and the on-shell scheme [18j. It is also worth noting that although 
two-loop corrections to the cross section may seem to be small, it is much harder to estimate 
their scale and behaviour for such a complicated observable as the parity-violating asymmetry 
to be measured by MOLLER experiment. 
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FIG. 1: Diagrams describing nonradiative M0ller scattering in the (1) t- and (2) u-channels. 



The two-loop EWC to the Born (~ M M +) cross section can be divided into two classes: the 
Q-part induced by quadratic one-loop amplitudes (~ MiM^), and the T-part corresponding 
to the interference of the Born and two-loop diagrams (~ 2ReMoM 2 l ^ oop ). The goal of this 
paper is to calculate the Q-part exactly. We show that the Q-part is much higher than the 
planned experimental uncertainty of MOLLER, which means that the two-loop EWC may be 
larger that previously thought. The large size of the Q-part demands a detailed and consistent 
consideration of the T-part, and that will be the next task of our group. 

II. GENERAL NOTATIONS AND MATRIX ELEMENTS 

Let us start by writing the cross section of polarized M0ller scattering with the Born kine- 
matics shown in Fig. [TJ 

e-(fci) + e"(pi) -> e"(fc 2 ) + e~(p 2 ), (1) 
such that, with the appropriate accuracy for the present paper, we find: 

3 3 

a = |~|M + Mx| 2 = ^-(M M + + 2ReM 1 M + + MiM^). (2) 

Here, a = da jd cos 9, 9 is the scattering angle of the detected electron with 4-momentum fc 2 
in the center-of-mass system of the initial electrons. The 4-momenta of initial (k\ and p\) and 
final (/c 2 and p 2 ) electrons generate a standard set of Mandelstam variables: 

8=(ki+ Pl ) 2 , t=(ki- k 2 )\ u={k 2 - Pl f. (3) 

It should also be noted that the electron mass m is disregarded wherever possible, in particular 
if m 2 s, —t, —u. 

Mo and M\ are the Born (0(a)) and one-loop (0(a 2 )) amplitudes (matrix elements), respec- 
tively. Let us describe the structure of M : 

M = M , t - M 0tU , M , u = M ,t(fc 2 p 2 ), M , t = V M/, Ml = i-IW jt J^, (4) 

z — / 7T" i" 
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where the t-channel upper and lower electron leg currents are 

JJ = u(/c 2 )7m(^ - a J 7sM^i), Jft = u^l^v 3 - a ] ^)u{pi). (5) 
The squared Born amplitude Mq forms the Born cross section: 

3 2 

a = — M M + = — [^ J (u 2 D u D jt + t 2 D iu D JU ) + X l fs 2 (D u + D iu )(D jt + D ju )]. (6) 
2s s 

A handy structure to use in the present study is 

D ir = — l — (i = >y,Z; r = t,u) 7 (7) 
r — mi 

which depends on the Z-boson mass mz or on the photon mass m 7 = A. The photon mass is 
set to zero everywhere with the exception of specially-indicated cases where the photon mass is 
taken to be an infinitesimal parameter that regularizes the infrared divergence (IRD). Another 
set of useful functions is 

\ i,k \ i,k \ i,k _i_ \ i,k \ i,k / c \ 

A± — A\ B Ai T ± A2 B A2 T i [o) 

These are combinations of coupling constants and Pbct), where Pb(t) are the degrees of polar- 
ization of electrons with 4-momentum k\ (pi). More specifically, 



1 B(T) - A V -PB(T)A A , M B (T) - A A -p B {T)A 



where 



\U = «V + aV, X l i = tfa? + aV, (9) 



v i = i, a 7 = o, v z = (J3 + 2s 2 w )/(2s w c w ), a z = I 3 J(2s w c w ). (10) 



The subscripts L and R on the cross sections correspond to Pb{t) = — 1 and Pb(t) = +1, where 
the first subscript indicates the degree of polarization for the 4-momentum k\ and the second 
one indicates the degree of polarization for the 4-momentum p\. Let us recall that I\ = —1/2 
and sw (c^y) is the sine (cosine) of the Weinberg angle expressed in terms of the Z- and W-boson 
masses according to the rules of the Standard Model: 



c w = m w /m z , s w = yl - & w . (11) 

We can present the one-loop amplitude M\ as a sum of boson self-energy (BSE), vertex (Ver) 
and box diagrams: 

Mi = M M - Mi,*, M ltU = M ljt (k 2 <-> p 2 ), 

M M = M BS E,t + M Ver , t + M Box , t . (12) 

We use the on-shell renormalization scheme from j2l|, 22|, so there are no contributions from 
the electron self-energies. The question of the dependence of EWC on renormalization schemes 
and renormalization conditions (within the same scheme) was addressed in our earlier paper 

Hi. 
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FIG. 2: One-loop t-channel diagrams for the M0ller process. The circles represent the contributions of self-energies and 
vertex functions. The u-channel diagrams are obtained via the interchange hi O P2- 



The infrared-finite BSE term can easily be expressed as: 

M BSEjt = i- ^ /;/;-'./"•'. (is) 

with 

£$' r = -D ir t^(r)D jr , (14) 

where Sy (r) is the transverse part of the renormalized photon, Z-boson and self-energies. 
The longitudinal parts of the boson self-energy make contributions that are proportional to 
m 2 /r; therefore they are very small and are not considered here. 

In order to get the electron vertex amplitude (2nd and 3rd diagrams in Fig. [2]), we use the 
form factors SF^ e A in the manner of paper [2lj], replacing the coupling constants v\ a? with 

form factors -> 5F^ Z)e , -> 5Fl {Z)e . Then, 

M Ver ,t = ( M ^>* + M J/H,t), M j/Bit = i^B^J^, M j/H , t = /^//,/>' / //"-'. (15) 

3=J>Z 

where the electron currents with vertices look like 

Bl = If t (v j -)• SF^, a? -»• 5F^ e ), = J W V -> <^y , -> (16) 

The infrared singularity is regularized by giving photon a small mass A and in the t-channel 
vertex amplitude can be extracted in the form: 

<r.* = ~ (log^-l)bg^Mo, t> (17) 

where e is the base of the natural logarithm. The rest (infrared-finite) part of t-channel vertex 
amplitude has the simple form The remaining (infrared-finite) part of the t-channel vertex 



amplitude has a simple form convenient for analysis and coding: 

,2 



M( ert = M Verii - = M Verii (A 2 -> s). 



The box term can be presented as a sum of all two-boson contributions: 



Mi 



Box,* 



M, 



77,t 



M lZ , t + M Z7)t + M zz , t + M W)t . 



We need to account for both direct and crossed 77, 7Z and ZZ-box.es: 

M ijtt = M° t + Mg t {,.., -./.). 

with Mg ft and Afj? t given by exact expressions in 4-dimensional integral form (4-point 
by 



(18) 

(19) 

(20) 
functions) 



M, 



D 

ij,t 



I - 



d 4 k 



MS, 



a 



(2tt) 2 J (k 2 -2k 1 k)(k 2 + 2p 1 k)((q-k) 2 -m 2 j )(k 2 -m 2 i ) 
xu(k 2 )^(v j - a J 7 5 )(fci - fc + m)Y(v l - a*7 5 )w(/ci) 
XttW7/i(f J - « J 75)(Pi + fc + m)^ v (v l - a l ^u{pi), 

d 4 k 



x 



(21) 



tt/ (2tt) 2 7 (k 2 - 2hk)(k 2 - 2p 2 k)((q - k) 2 - m 2 )(k 2 - m 2 ) 

xu(k 2 )'j fJ '(v j - a j j 5 )(ki - fc + m)^ v (v % - a % ^^)u{ki) x 
X«(p 2 )7^K - « l 75)(P2 -k + m^^v 3 - a>^ b )u{pi). 



(22) 



Obviously, for PUPU-boxes we only need the crossed expression ( 1221) . 
The infrared parts of the 77- and 7Z-boxes in the t-channel are similarly given by 



^77( 7 Z+Z 7 ),i 



a/1 -u 

- U lo s — lo S 

7T \2 S 



— U , —US 7T S\ 

_ + _ + !7r l 0g _)M l -'< ». 



(23) 



Using asymptotic methods, we can significantly simplify the box amplitudes containing at least 
one heavy boson (see, for example, [14J, where simplifications were done on the cross-section 
level). Then 



M 



/ax 2 
) X 



M* Ztt + M f z ^ t = [M lZ>t + M Zl ^ - (m 7 \, + 
2 

+ log ^y^u(k 2 )^(v z - a Z 7 5 )(-7 a bX/i;i) • u{p 2 )l»{v z - a z lh ) loilv u{j>x) + 



3 
2 



+ 



2 

+ log — Wfc 2 )7 M (^ - a z 75 ) 7 Q 7^(A; 1 ) ■ u{p 2 ) lulalll {v z - a z l5 )u( Pl ) 
— u / 



(24) 



M 



-a\ 2 1 
-il - 



W 16m 2 z 



u{k 2 )Y(v B - a B 7 5 )(-7 a )7 I/ M(A; 1 ) • u(p 2 )^{v B - a B l5 ) lalv u{p x ) + 



+ M(fc 2 ) 7 ^(^ - a B 75 )7 a 7^(fc 1 ) • u{p 2 ) lvlalll {v B - a B l5 )u( Pl ) 



B „B, 



(25) 
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Mi 



ww,t — —i I 



/a\ 2 1 



\7r/ 16m^ 



u(k 2 )j' x (v c - a c >y 5 )>y a Yu(k 1 ) ■ u(p 2 hulal^(v u - a u j 5 )u(p 



c „c. 



with the coupling-constants combinations for ZZ- and VTiy-boxes 



^ Z ) 2 + (a 



Zn 2 „B 



2w z a z , = a°' = l/(4s^) 



, (26) 



(27) 



Now we are ready to present the one-loop complex amplitude as the sum of IR and IR-finite 
parts: 



Afi = Mi + M{, Mi = -h^Mo, M{ = Mbse + M( ei + M| ox + M a , 

TV 2 



where 



** = 4Blog-L 



and the complex value B can be presented in form (see, for example, 23] ) 

tu 



B = lot 



m 2 s 



1 + Z7T. 



The amplitudes from the non-factorized part of the boxes are given by 



ReM n = - 



a 
2^ 



(L 2 U + 7r 2 )M 0it - [L\ + n 2 )M 0tU 



where L r . = log(— s/r). 



III. EXTRACTION OF INFRARED DIVERGENCES 



(2* 



(29) 



(30) 



(31) 



Now we should make sure that the infrared divergences are cancelled. In a similar way as 
it was done for amplitudes, we present the complex interference term ai and differential cross 
section gq as sums of A-dependent (IRD-terms) and A-independent (infrared-finite) parts: 



<7i = —MiM+ = <7* + a{, o Q = ^—MiMt = o X Q + cr f Q . 



(32) 



The one-loop cross section which we denote <7i = Re<7i was carefully evaluated with full control 
of the uncertainties in paper [lij • The term ctq (see (j2J)) is called the Q-part of the two-loop 
EWC and is the main subject of the present paper. 

If we substitute the amplitudes derived in Section II to the left-hand- side of (j2J), and compare 
the result with the right-hand side of this equation, we will get the same expression for o\ as 



given in [14J. The simplest form for o\ (see formula (42) of |14|) is then: 



a 



A « cA^ 
7T 



(33) 



The infrared-finite part erf can be conveniently to presented via the relative dimensionless 
correction: 



a[ = -S(a . 



71 



(34) 



After some transformations, the value crA is given by 



4 \7T/ 



(35) 



Finally, the infrared-finite part Oq expressed via the relative dimensionless corrections has form 



7T 

2s 



M{M{ + =Q 2 6 f Q a . 



(36) 



IV. BREMSSTRAHLUNG AND CANCELLATION OF INFRARED DIVERGENCES 



To evaluate the cross section induced by the emission of one soft photon with energy less 
then u, we follow the methods of [Hf (see also 25]). Then this cross section can be expressed 
as: 



cy-=aj + al (37) 

where a\ 2 have the similar factorized structure based on the factorization of soft-photon 
bremsstrahlung: 

°] 

7T 



a J = -Re [-5* + RAao, a\ = -Re[(-# + i?i)*<Ti] ; 



(3f 



where 



i2i = -45 log 



log 2 



7T 



log' 



(39) 



2w u em 2 3 t 

The first part of the soft-photon cross section, aj, cancels the IRD at the one- loop order, 
while the second part, aj, cancels the IRD a the two- loop order, with half of al going to the 
cancellation of IRD in the Q-part and the other half going to treat IRD in the T-part: 

„ 1 



a 



T 



-OA 



(40) 



To obtain the term —b\ + R\ in Eq. 
the phase space of one real soft photon. It can be done according to 24 

i r d 3 k 



we must calculate the 3-dimensional integral over 

in c.m.s: 



-Sl + Ri =L(\,u) 



4tt J k 



■T?(k)Tp(k), 



where 



h a b a n a n a 

T a (k) - 



(41) 



(42) 



k\k k 2 k pik p 2 k 

The difference between the estimation relying on the soft part only and the result obtained by 
separation into the soft and hard parts at lowest order is rather small (see 14J), so we believe 
that the soft cross section will provide the sufficient accuracy at second order as well. 

At last, the cross section induced by the emission of two soft photons with a total energy less 
then u can be written as: 

2 



a 



77 



a 
2 \n 



R2 O"0 



(43) 
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where | is a statistical factor caused by the indistinguishability of two final photons and R2 = 
^7i 2 \B\ 2 . The detailed calculations of cr 77 are shown in Appendix A. 

Just like a 7 , the cross section a 77 is divided into equal halves, with a half going to cancel the 
IRD in the Q-part and a half going to the T-part: 

aQ = = \a^. (44) 

Combining all the terms together, we get the infrared-finite result at both the first and second 
orders. The first (second) order is given by the first (second) term on the LHS of the equation 
below: 

Rela.+a^ + iaQ + al + a^) = ^R e [R 1 + 5[]a + (^Re^tf + ^ + 1^- ±R 2 ]a . (45) 



7T 



NUMERICAL RESULTS 



For the numerical calculations we use a = 1/137.035999, m w = 80.398 GeV, and m z = 
91.1876 GeV as input parameters in accordance with j26[. The electron, muon, and r-lepton 



masses are taken to be m e = 0.510998910 MeV, m M = 0.105658367 GeV, m T = 1.77684 GeV, 
while the quark masses for vector boson self-energy loop contributions are taken to be m u = 
0.06983 GeV, m c = 1.2 GeV, m t = 174 GeV, m d = 0.06984 GeV, m s = 0.15 GeV, and m b = 
4.6 GeV. The values of the light quark masses were extracted using the fact that they provide 
shifts in the fine structure constant due to hadronic vacuum polarization Aa^(m|)=0. 02757 



27 , where 



f—„. A v V / / 



f=u,d,s 



and Qf is the electric charge of fermion / in proton charge units q (q — y/ Ana) . 

On the other hand, the contribution of hadronic vacuum polarization to the fine structure 
constant also can be evaluated using the dispersion relation: 

00 

where V means that the principle value of the integral should be considered and cr/i(s) is the 
cross section of hadron production in e + e~ annihilation. In the case of small energies this cross 
section can be approximated by the cross section of the pion production channel e + e~ —> 7r + 7r~: 



*M = f^, A, = V 1 - — • < 48 > 

3s V s 

thus giving the following contribution to Aa[ 5 ^(s): 

Using Eq. (I46p and Eq. (1491) we can incorporate the use of the light quark masses as parameters 
regulated by the hadronic vacuum polarization in our calculations. 



0.03 




FIG. 3: Virtual and bremstrahlung contributions to the relative correction to unpolarized cross section vs. the photon 
mass A at 9 = 90° and Si ab = 11 GeV. 



Finally, for the mass of the Higgs boson, we take tuh = 115 GeV. Although this mass is 
still to be determined experimentally, the dependence of EWC on m# is rather weak. For the 



maximum soft photon energy we use u = 0.05y/s, according to 14J and 

Let us define the relative corrections to the Born cross section due to a specific type of 
contributions (labeled by C) as 

5 C = (a c -a°)/a°, C = 1-loop, Q, T, ... 

The parity-violating asymmetry is defined in a traditional way, 

. v LL + <J LR — a RL — cr RR Cll — Crr / ca x 
Alr = = , (50) 

<?LL + <?LR + ORL + ORR & LL + ^ LR + &RR 

and the relative correction to the Born asymmetry due to C-contribution is defined as 

£>A = i A LR ~~ A° LR )/A° LR . 

Fig. El plotted for 9 = 90° and = 11 GeV, clearly demonstrates that the relative correction 
to the unpolarized cross section is independent of the photon mass A. We can also see a 
quadratic dependence on the log scale of A for both the virtual and bremstrahlung contributions. 

The left frame of Fig. |4] depicts the relative corrections to the asymmetry at = 11 GeV 
versus the scattering angle 8 in c.m.s. The lower line shows the corrections to the asymmetry 
with only one-loop EWC taken into account, and the upper line shows the combined one-loop 
and Q-part corrections. As expected, both of them are symmetric along the line 9 = ir/2, have 
a minimum at 9 = 90°, and depend on the scattering angle quite weakly. 

The difference of these two effects is an absolute correction defined by 



A — t yjl-loop+Q ^0 W/lO (a 1-loop a0 \ / a0 _/ zl l-loop+Q ^l-loopw^O 

— {S±LR ~ AlRII^LR ~ \ A LR ~ A LR> I A LR ~ \ A LR ~ LR ) I A L 



LR 
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FIG. 4: The relative corrections to the asymmetry (left) and the absolute correction Aa (right) vs scattering angle 9. 



and depicted in the right frame of Fig. HI Here we can see that the Q-part gives quite a 
significant contribution, with Aa reaching a maximum of 0.0419 at 9 = 90°. Taking into 
account that MOLLER's planned experimental error to the PV asymmetry is ~ 2% or less, we 
see that it is necessary to continue to work on the two-loop EWC, staring from the T-part. 

Fig. [5] shows the relative (labeled as 1-loop and 1-loop+Q) corrections and absolute 
corrections (labeled by Q) versus yfs at 6 = 90°. In the high-energy region (^/s > 50 GeV) our 



one- loop result (see [14J) is in excellent agreement with the result from [28[ if we use the same set 
of Standard Model parameters. As one can see from Fig. |5l the scale of the Q-part contribution 
in the low-energy region is approximately constant, but grows sharply at y/s > rriz, where the 
weak contribution becomes comparable to the electromagnetic. This increasing importance of 
the two-loop contribution at higher energies may have a significant effect on the asymmetry 
measured at future e~e"-colliders. 



VI. CONCLUSIONS 



Experimental investigation of M0ller scattering is not only one of the oldest tools of modern 
physics, but also a powerful probe of new physics effects. The new ultra-precise measurement 
of the weak mixing angle via 11 GeV M0ller scattering planned at JLab (MOLLER) - as 
well as experiments proposed at future high-energy electron colliders - will require that the 
higher-order effects to be taken into account with the highest precision possible. 

In this work, we build on the study of the one-loop electroweak radiative corrections to the 
cross-section asymmetry of the polarized M0ller scattering at 11 GeV initiated by our group in 



14] . and address some of the two-loop effects. At this stage, we perform a detailed calculation 
for the part of the two-loop electroweak radiative correction induced by squaring one-loop 
diagrams. 
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FIG. 5: Relative (labeled by 1-loop+Q and 1-loop) and absolute (labeled by Q) corrections to PV-asymmetry vs i/s. The 
filled circle corresponds to our predictions for the MOLLER experiment. 



The two-loop EWC to the Born (~ MqMq) cross section is divided into the T-part, which is 
the interference of Born and two-loop diagrams (~ 2ReM M 2 l l ioop ), and the Q-part, induced by 
quadratic one-loop amplitudes (~ M\M^), which we evaluate here. The results are presented in 
both numerical and analytical form, with the infrared divergence explicitly cancelled. Also, we 
clearly demonstrate the important role of the imaginary part of amplitude, which is consistently 
taken into consideration both in the infrared-finite and divergent terms. 

As one can see from our numerical data, at the MOLLER kinematic conditions, the part of 
the NNLO EWC we considered in this work can increase the asymmetry by up to ~ 4%. The 
corrections depend quite significantly on the energy and scattering angles; at the high-energy 
region of y/s ~ 1000 GeV achievable in the planned experimental program of the ILC, the 
estimated contribution of the quadratic EWC can reach +14%; for 3 TeV at CLIC, it would be 
+42%. We see that the large size of the Q-part demands detailed and consistent consideration 
of the T-part, which will be the next task of our group. It is impossible to say at this time 
if the Q-part will be partially enhanced or cancelled by other two-loop radiative corrections, 
although it seems probable that the two-loop EWC may be larger than previously thought. 
Although an argument can be made that the two-loop corrections are suppressed by a factor of 
air relative to the one-loop corrections (see [HI], for example), we are reluctant to conclude that 
they can be dismissed, especially in the light of 2% uncertainty to the asymmetry promised by 
MOLLER. 

Since the problem of EWC for the M0ller scattering asymmetry is rather involved, a tuned 
step-by-step comparison between different calculation approaches is essential. One of the im- 
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portant results of this work is the correctness of our calculations, which was controlled by a 
comparison of the results obtained from the equations derived by hand with the numerical data 
obtained by a semi-automatic approach based on FeynArts, FormCalc, LoopTools and Form. 



These base languages have already been successfully employed in similar projects ([14J, |19j). 
so we are highly confident in their reliability. 

In the future, we plan to address the remaining two-loop electroweak corrections which may be 
required by the promised experimental precision of the MOLLER experiment and experiments 
planned at ILC. 
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FIG. 6: Double-bremsstrahlung diagrams for M0ller scattering in t-channel corresponding to M\ x . The u-channel diagrams 
are obtained by interchange k-2 o P2 ■ 



Appendix A: Detailes of calculations for the case of emission of two real soft photons 



First, let us calculate the amplitudes corresponding to the emission of two real soft photons 
(see Fig. Ej), 

e~(A;i,0 + e-(pi,77) ->e~(A; 2 )+e-(p 2 )+7(A;)+7(p) (Al) 

in t- and u-channels with z-boson exchange (i = 7, Z). The amplitudes are labeled as 
M-q, M{ 2 i Mis, ■■■■> where the first (second) subscript denotes the origin of the first j(k) (second 
j(p)) emitted photon: 1 - emitted from electron e~(ki), 2 - from electron e~(k 2 ), 3 - from 
electron e~(pi) and 4 - from electron e~(p 2 ). The exact expression for M\ Y is the following: 

M* u = i{2nefN kl N k2 N Pl N P2 N p N k ■ u(k 2 )Y(v l - a t l5 ) r ^ 7 a e a (p) r ^—^e p (k)u(k l ) ■ 

fc 4 — m k 3 — m 

v{p2)lu{v l - a l j 5 )v(p!) ■ — 2 ■$(k 1 +p 1 -k 2 -p 2 -k-p), (A2) 
where N k = ( 27r )3/2 Using the Dirac equation and taking k — > 0, we can simplify 

1 kx-k + m B fci + m « . . 

-Y u {ki) = J-, pr^ ^ Pu \ki) ~ — — — -7 u(fci) = 



k 3 -m (h - k) 2 - m 2 -2k x k 

= Z^( 2fc f + Vhfe + = "^M- ( A3 ) 



Analogously, at fc,p — > 0, 



-iMh) = , „ k \ , , u(fci). (A4) 



A; 4 -m -fci(fc + p) + fcp 



Finally, the amplitude MJ X at — >■ has the following form, with the convenient factoriza- 
tion from the Born amplitude: 

k a kP 

^'■ £ -«*'H*ira' M; (A5) 
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In the same manner, we get 

M^kp^o = e 2 N p N k ■ e a (p)ep(k) 

M[ 2 \ k)P ^ = e 2 N p N k ■ e a (p)ep(k) 

Mi 3 \ k , p ^o= e 2 N p N k - e a (p)ep(k) 

u\k, P ^o = e 2 N p N k ■ e a {p)ep(k) 

3 4lfc,p^o = e 2 N p N k ■ e a (p)ep(k) 

Ml 3 \ k>p ^ = e 2 N p N k ■ e a (p)ep(k) 

MJ 4 |/;^o = e 2 N p N k ■ e a (p)ep(k) 

M 23\k, P ^o = e 2 N p N k ■ e a {p)e p {k) 

M l 24 \ k:P ^o = e 2 N p N k ■ e a (p)ep(k) 



h a h p 



M, 



m: 



k 2 p)(k 2 (k + p) + kp) °' 

<^2 "-1 



- • M l 
k 2P )(-k 1 k) °' 



a P 

ViVi 



-pik)(-p 1 (k + p) + kp) 

M, 



• M, 



o- 



a P 
P2P2 



P2p)(p2(k + p) + kp) 
■Mi 



o- 



a P 
P2P1 



P2P){-Pik) 
k*p{ 



-k lp )(- pl k) ' M °' 

-kip)(p 2 k) 
k a 2 p{ 



hp)(-Pik) M °' 



■ M l 

k 2 p)(P2k) °' 



(A6) 



Now we need to sum the terms generated by the substitution k p. For the 11-, 22-, 33-, 
44-cases it works as the following: 



M[ 1 \ Kp ^ + (k^p) = e 2 NpN k M^— 

= e 2 N p N k M l 



e a (p)e /3 (k)kfk f 1 



+ 



e a (k)ep(p)k?k^ 



k)(—kx (k + p) + kp) (-kip)(-k x (k + p) + kp) 
1 1 \ e a (p)ep(k)k?k? 



e N p N k M, 



k\k —kip/ (—ki(k + p) + kp) 

i e a (p)e [S (k)k^k^ 



(hk)(k lP ) ■ 

As a result, the total i(w)-channel amplitude is given by 

M t y fe ,^o = e 2 N p N k ■ e a (p)ep(k) ■ T a (p)T?(k) ■ M^ {u) . 
Then, the cross section of (1A1I) has the form 

d s kd 3 p ■ (eN k ) 2 (eN p ) 2 T a (p)T a (p)T^(k)Tp(k) 
d 3 k d 3 p 



(A7) 



(A8) 



a 
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a °2 



a\ 2 1 / 1 \ 2 



7V J 2V47T/ 



fco+po<" 



ko Po 



T a (p)T a (p)T?(k)Tp(k). 



(A9) 
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It is possible to prove that 

1 \2 f d 3 kd 3 p 



(i-) 2 j ^ d ^-T a (p)T a (p)T^(k)T,(k) = \L(X,u)\ 2 -R 2 . (A10) 

k +po<u 

If we change the condition ko + po < u to simply k < u, p < u, the term R 2 will go to zero. 
However, the conditions ko < u, po < u are not valid from the experiment point of view. 
Let us calculate R 2 exactly. First, we introduce the notation 

I(a 1 ,a 2 ;b 1 ,b 2 )=(±-) 2 [ • T*(p)T a (p)T?(k)T p (k). (All) 

\AixJ J k p 

a 1 <k <a,2, b 1 <p <b 2 

From Eq.f l4ip it is obvious that 

7(0, a; 0,6) = L(X,a)L(X,b)* , (A12) 

which, of cause, works at A <C a, b only. At da, db <^ a, b, the simple geometry considerations 
give the equation 

I (a, a + da;b,b + db) = 1(0, a + da;0,b + db) + 1(0, a; 0, b) - 1(0, a;0,b + db) - 1(0, a + da; 0, b). 

Simplifying using Eq.( 1A12| ) and Eq.( )4~Tj) . we get 

ti 77, 77\ lHn i2, a + da . b + db ^.^,-dadb ,. ... 

I(a, a + da;b,b + db) = 1Q\B\ 2 log log — - — « I6|5| 2 -. (Al3) 

a b a b 

Finally, comparing Eq. flAlOp and Eq. flAl3p . we conclude: 

R 2 = Yl(a,a + da;b,b + db) = lG\B\ 2 [— [ — = l6|5| 2 Li 2 (- )" = -tt 2 \B\ 2 . (Al4) 
n J a J b \uJ o 3 

il u-a 



Here, f2 = {a < cu} R {b < lu} fl {a + b > u}. Our result for R 2 exactly agrees with [23 



